Adiabatic loading of bosons into optical lattices 
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The entropy-temperature curves are calculated for non-interacting bosons in a 3D optical lattice and a 2D 
lattice with transverse harmonic confinement for ranges of depths and filling factors relevant to current experi- 
ments. We demonstrate regimes where the atomic sample can be significantly heated or cooled by adiabatically 
changing the lattice depth. We indicate the critical points for condensation in the presence of a lattice and show 
that the system can be reversibly condensed by changing the lattice depth. We discuss the effects of interactions 
on our results and consider non-adiabatic processes. 
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I. INTRODUCTION 

Neutral bosonic atoms in optical lattices have been used 
to demonstrate quantum matter- wave engineering QJ 0], 
the Mott-insulator quantum-phase transition |3], and explore 
quantum entanglement [4]. The many favorable attributes of 
optical lattices, such as the low noise level and high degree of 
experimental control, make them an ideal system for imple- 
menting quantum logic Ollg]. A central part of many of these 
proposals is the use of the Mott-insulator transition to prepare 
the system into fiducial state with precisely one atom per site, 
and negligible quantum or thermal fluctuations. 

The usual path for preparing a sample of quantum degener- 
ate bosons in an optical lattice consists of first forming a cold 
Bose-Einstein condensate in a weak magnetic trap, to which a 
1, 2 or 3D lattice potential is adiabatically applied by slowly 
ramping up the light field intensity. Ideally this process will 
transfer a condensate at zero temperature into its many-body 
ground state in the lattice. Of course condensates cannot be 
produced at T = OK, yet to date the role of temperature has re- 
ceived little attention even though it may play a crucial role in 
the many-body properties of these systems. Indeed, due to the 
massive energy spectrum changes the system undergoes in lat- 
tice loading, the initial and final temperatures are not trivially 
related. It is also of interest to understand how T c changes 
in the lattice to assess the effect the periodic potential has on 
condensation. Relevant to these considerations, experiments 
reported in 1 7] examined evaporative cooling of atoms in a 
combined magnetic trap and ID optical lattice, and showed a 
significant decrease in the critical temperature for a relatively 
shallow lattice depth. 

Using entropy comparison Olshanni and Weiss |8] have 
considered how a thermalized system of bosons in an optical 
lattice would be transformed through adiabatic unloading into 
simple traps, with a view to producing a condensate optically. 
Their approach takes into account spatially inhomogeneous 
potentials superimposed upon the lattice, but they assume that 
each lattice site is occupied by no more that one boson and 
that the tunneling rate between sites is zero — strictly valid 
only for infinitely deep lattices. 

In this paper we begin by considering the thermodynamic 
properties of an ideal gas of bosons in a 3D cubic lattice. 
Working with the grand canonical ensemble we use the ex- 
act single particle eigenstates of the lattice to determine the 



entropy-temperature curves for the system for various lattice 
depths and filling factors. We use these curves for analyz- 
ing the effect of the loading process on the temperature of the 
system, and show that sufficiently cold atomic samples can be 
significantly cooled through loading into a lattice. By analyz- 
ing the nature of the energy spectrum we explain this counter- 
intuitive notion that adiabatic compression of a system can 
lead to cooling. 

Interactions between particles are not accounted for in our 
calculations of the thermodynamic properties, however by 
considering the Bogoliubov excitation spectrum in the lattice 
we discuss the modifications interactions should introduce to 
our results. We also address how robust our predictions are to 
non-adiabatic effects in the loading procedure. 

In the last part of this paper we investigate the thermody- 
namic properties of bosonic atoms in a two-dimensional lat- 
tice. An important experimental consideration in this case 
is that the intensity envelope of the lasers used to make the 
optical lattice gives rise to an additional slowly-varying po- 
tential perpendicular to the plane of the lattice sites. For the 
case of a red detuned lattice this potential is confining and 
approximately harmonic in the region where the atoms are 
trapped. With increasing laser intensity the lattice (ground 
band) degrees of freedom become more degenerate, whereas 
the energy spacing of harmonic degrees of freedom increases 
due to the strengthened transverse confinement. Competition 
between these two effects considerably modifies the nature 
of the heating or cooling that occurs during lattice loading. 
We model the thermodynamic properties of this system and 
present numerical calculations for the entropy-temperature 
curves for parameter regimes relevant to current experiments. 



II. FORMALISM - 3D LATTICE 

A. Single Particle Eigenstates 

We consider a cubic 3D optical lattice made from 3 in- 
dependent (i.e. non-interfering) sets of counter-propagating 
laser fields of wavelength A, giving rise to a potential of the 
form 
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where k — 2tt/X is the single photon wavevector, and V is the 
lattice depth. We take the lattice to be of finite extent with a 
total of N s sites, consisting of an equal number of sites along 
each of the spatial directions with periodic boundary condi- 
tions. The single particle energies e q are determined by solv- 
ing the Schrodinger equation 
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for the Bloch states, ipq(r), of the lattice. For notational sim- 
plicity we choose to work in the extended zone scheme where 
q specifies both the quasimomentum and band index of the 
state under consideration. By using the single photon recoil 
energy, Er = h 2 k 2 /2m, as our unit of energy, the energy 
states of the system are completely specified by the lattice 
depth V and the number of lattice sites N s (i.e. in recoil units 
e q is independent of k). 

It is useful to review the tight-binding description of the 
ground band. Valid for moderately deep lattices, the tight- 
binding model only accounts for nearest neighbor tunneling, 
and leads to the analytic dispersion relation 



TB 



3 — cos(qjd) 



(3) 



where m* = f h 2 V q e q 
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q=0 



is the effective mass at 



q = 0, and a = A/2 is the lattice site spacing. 
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Figure 1: Density of states for a N s ~ 3 x 10 4 site cubic lattice at 
various depths. For a depth of approximately V = 2Er a gap devel- 
ops in the density of states. In (c) we illustrate the energy gap, e gap , 
and ground bandwidth, ebw- Points are determined by numerically 
averaging the exact spectrum over a small energy range. 



where fi is found by ensuring particle conservation. The en- 
tropy of the system can then be expressed as 

S = k B (log Z + 0E - f i(3N p ) , (5) 

where (3 — 1/ k B T, and E = — d In Z /d(3 is the mean energy. 



B. Equilibrium Properties 

Our interest lies in understanding the process of adiabati- 
cally loading system of N p bosons into a lattice. (The require- 
ments for adiabaticity in this system are not well understood, 
though the timescales for adiabatically changing the lattice 
depth are expected to become long in deep lattices, which we 
discuss further in Sec. IHI E> . Under the assumption of adia- 
baticity the entropy remains constant throughout this process 
and the most useful information can be obtained from know- 
ing how the entropy depends on the other parameters of the 
system. In the thermodynamic limit, where N s — > oo and 
N p — > oo while the filling factor n = N p /N s remains con- 
stant, the entropy per particle is completely specified by the 
intensive parameters T, V, n. The calculations we present in 
this paper are for finite size systems, that are sufficiently large 
to approximate the thermodynamic limit. 

To determine the entropy, the single particle spectrum {e q } 
of the lattice is calculated for given values of N s and V. We 
then determine the thermodynamic properties of the lattice 
with N p bosons in the grand canonical ensemble, for which 
we calculate the partition function Z 

logZ = - log (l - e-^o"")) , (4) 
q 



1. Density of states 

How the thermodynamic properties of a system of bosons 
change as they are adiabatically loaded into a lattice intimately 
reflects how the lattice modifies the microscopic energy spec- 
trum. In this regard the density of states function affords con- 
siderable insight into the behaviour of the system. In Fig. ^ 
we illustrate how the distribution of available single particle 
states changes for various lattice depths. 

In general we note that the lattice leads to a substantial 
change in the the density of states, g(e), for the system. In 
the absence of the lattice (Fig. [TJ a )) the density of states is 
that for free particles in a box, and is proportional to y/e. The 
smoothness of g(e) is disrupted by the presence of a lattice, 
which causes flat regions in the energy bands giving rise to 
peaked features in the density of states, known as van Hove 
singularities \9]. The van Hove singularities in the first and 
second energy bands are clearly visible in Figs, ^b)-(d). For 
sufficiently deep lattices an energy gap, e gap , will separate the 
ground and first excited bands. For the cubic lattice we con- 
sider here, a finite gap appears at a lattice depth of V ~ 2Er 
fllil (see Fig. ^b)), and beyond this depth the gap increases 
with lattice depth (see Figs. HcMd))- In forming the gap 
higher energy bands are shifted up-wards in energy, and the 
ground band becomes compressed — a feature characteristic 
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Figure 2: Entropy versus temperature curves for a N s ~ 3 x 10 4 site 
cubic lattice, with filling factor n = 0.25 at various depths V = 
to 20Er (with a spacing of 2Er between each curve). The entropy 
plateau So is shown as a dashed line and the condensation point is 
marked on each curve as a circle. Dotted line marked A shows a path 
along which adiabatic loading into the lattice causes the temperature 
to increase. Dotted line marked B shows a path along which adia- 
batic loading into the lattice causes the temperature to decrease. 



of the reduced tunneling between lattice sites. We refer to 
the energy range over which the ground band extends as the 
(ground) band width, £bw- This quantity decreases in mag- 
nitude exponentially with V (see Figs^b)-(c)), causing the 
ground band to have an extremely high density of states for 
deep lattices; generally in this limit quantum many-body ef- 
fects will significantly modify the energy spectra of the lowest 
band from the non-interacting states we use here. We will dis- 
cuss the effects of interactions in Sec. IIII Dl 



III. RESULTS 

In Figs. I2l4l we show entropy-temperature curves for vari- 
ous lattice depths and filling factors n. These curves have been 
calculated for a lattice with 31 lattice sites along each spatial 
dimension, i.e. N s w 3 x 10 4 . The condensation tempera- 
ture is defined as that at which 1 % of all particle occupy the 
ground state, and is indicated in Figs. 12141 We note that being 
a finite system the transition is not discontinuous, however the 
onset of condensation is rapid and changing the requirement 
to 5% makes little observable difference in the critical point 
locations. 

An important common feature to these curves is the dis- 
tinct separation of regions where adiabatic loading causes the 
temperature of the sample to increase or decrease, which we 
will refer to as the regions of heating and cooling respectively. 
These two regions are separated by a common point that all 
curves approximately pass through, and we denote by its co- 
ordinates as {Tq, Sq}. The reason for the existence of this 
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Figure 3: Entropy versus temperature curves for a N B ~ 3 x 10 4 
site cubic lattice, with filling factor n = 1 at various depths V = 
to 2QEr (with a spacing of 2Er between each curve). The entropy 
plateau So is shown as a dashed line and the condensation point is 
marked on each curve as a circle. 
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Figure 4: Entropy versus temperature curves for a N s ~ 3 x 10 4 
site cubic lattice, with filling factor n = 4 at various depths V = 
to 20Er (with a spacing of 2Er between each curve). The entropy 
plateau So is shown as a dashed line and the condensation point is 
marked on each curve as a circle. Dotted line marked C shows a path 
along which an initially Bose-condensed system uncondenses as it is 
loaded into the lattice at V « 4Er. 



point will be discussed below. For cases considered in Figs. 
I2I4I Tn is in the range 0.5 — 1-Er, and to clearly indicate the 
vertical separation of the heating and cooling regions, we have 
marked Sq as a horizontal dashed line. 

We now explicitly demonstrate the temperature changes 
that occur during adiabatic loading using two possible adia- 



4 



batic processes labeled A and B, and marked as dotted lines 
in Fig. |2 Process A begins with a gas of free particles in 
an initial state with S > So, T > Tq. As the gas is loaded 
into the lattice the process line indicates that the temperature 
increases rapidly with the lattice depth. Conversely process 
B begins with a gas of free particles in an initial state with 
S < Sq, T < Tq, and the lattice loading causes a rapid de- 
crease in temperature. This behavior can be qualitatively un- 
derstood in terms of the modifications the lattice makes to the 
energy states of the system. As is apparent in Fig. ^ the 
ground band flattens for increasing lattice depth causing the 
density of states to be more densely compressed at lower en- 
ergies. Thus in the lattice all these states can be occupied at a 
much lower temperature than for the free particle case. As we 
discuss below, So is the maximum entropy available from only 
accessing states of the lowest band, and so for S > So higher 
bands are important. As the lattice depth and hence e gap in- 
creases, the temperature must increase for these excited states 
to remain accessible. 



A. Entropy plateau 



1. ksT <g; EBW- 

When the temperature is small compared to the ground 
band width, only energy states near the ground state of the 
lowest band are accessible to the system. These states exhibit 
a quadratic dependence on the magnitude of the quasimomen- 
tum and the entropy of this system is well described by the 
free particle expression if the bare particle mass is replaced 
by the effective mass i.e. m — ► m* . In general this regime oc- 
curs only at very low temperatures, and except for extremely 
low filling factors the system will be condensed in this region. 
In this limit the expression for entropy corresponds to the ex- 
pression for a condensed gas of free particles with mass m* 



S 



-AU 3 C(5/2) 



n*kBT 
2irh 2 



(8) 



where £(n) = J2T=i V^" I s tne Riemann Zeta function and 
N s a 3 is the volume of the system fioll . We note that in 
this regime the critical temperature for condensation scales as 
T c ~ (m*) -3 / 2 , so that T/T c is independent of m* and hence 
the lattice depth. 



In Figs. 12141 a horizontal plateau (at the level marked by the 
dashed lines) is common to the entropy-temperature curves 
for larger lattice depths (V > 8Er). This occurs because for 
these lattices eew "C £gap, and there is a large temperature 
range over which states in the excited bands are unaccessible, 
yet all the ground band states are uniformly occupied. The 
entropy value indicated by the dashed line in Figs. l2l4l corre- 
sponds to the total number of A^-particle states in the ground 
band. Since the number of single particle energy states in the 
ground band is equal to the number of lattice sites, the to- 
tal number of available Ap-particle states, which we define as 
Qq, is the number of distinct ways N p identical bosons can be 
placed into N s states, i.e. 



O = 



(n p + n s -i)\ 



N P \(N S - 1)! 



(6) 



The associated entropy So 



approximation 



Sq 



N p \og 



N p + N s -1 



(7) 



As described earlier, this entropy value separates the heating 
and cooling regions for the lattice. 



B. Scaling 

Here we give limiting results for the entropy-temperature 
curves. 



2. Tight-binding limit with feT< e gap 



As discussed in Sec. IIII Al when the temperature is small 
compared to the energy gap only states within the ground 
band are accessible to the system. In addition when the tight- 
binding description is applicable for the initial and final states 
of an adiabatic process, the initial and final properties are re- 
lated by a scaling transformation. 

To illustrate we consider our initial system to be in equilib- 
rium with entropy Si < So, in lattice of depth Vi which we 
take to be sufficiently deep enough for tight-binding expres- 
sion Q to be a good description of the ground band energy 
states. If an adiabatic process is used to take the system to 
some final state at lattice depth V) (also in the tight-binding 
regime) it is easily shown that the macroscopic parameters of 
the initial and final states are related as 



Sterling's 




= aEi, 


(9) 






= aTi, 


(10) 




f'f 




(11) 



where the scale factor 



m*/m 



/ 



(12) 



is the ratio of the effective masses, and the subscripts i and 
/ refer to quantities associated with the initial and final states 
respectively. This type of scaling relation leaves the occupa- 
tions of the single particle levels (including the condensate 
occupation) unchanged. This means that being adiabatic does 
not require redistribution through collisions and may allow the 
lattice depth to be changed more rapidly. However, adiabatic - 
ity with respect to the many-body wavefunction will become 
a significant constraint as the lattice depth increases and will 
become the limiting timescale. 
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C. Condensation 

The temperature and entropy at which bosons condense 
generally changes with lattice depth as indicated in Figs. 12141 
We note that for high filling factors the condensation points for 
different lattice depths occur over a wide range of entropies, 
suggesting that the degree of condensation will be greatly af- 
fected by adiabatic lattice loading. 

For instance, consider the adiabatic process indicated by 
the dotted line and labeled C in Fig. |4] The system starts 
as a Bose-condensed gas of free particles. However, as the 
lattice depth increases the condensate fraction decreases until 
the system passes through the transition point and becomes 
uncondensed. This process is reversible, and is analogous to 
the experiments by Stamper- Kurn et al. lHHl . where a Bose 
gas was reversibly condensed by changing the shape of the 
trapping potential. 

We note that the free particle case with n = 1 shown in 
Fig. [3] condenses at an entropy approximately equal to the 
plateau entropy So, given by Eq. Q. For an ideal gas of free 
particles condensation occurs when the entropy of the sys- 
tem is S cond = ffcu-WpCCD/CCf )> (e.g. see OjJ|), whereas 
the plateau entropy for a lattice with unit filling is So(n — 
1) = 2k B N p \og{2), (obtained by setting (N s - 1) -> N s 
and N p = N s in Eq. (0). In comparing the values we find 
So(n = lVScond ~ 1-0798, in agreement with numerical 
observation. For fixed N s the filling factor is proportional 
to N p , thus the free particle condensation entropy scales as 
Scond ~ n whereas the entropy plateau goes as So ~ log(n). 
So we conclude that for n > 1 condensation occurs (for free 
particles) above the plateau, whereas for n < 1 it occurs be- 
low the plateau. 

D. Interaction effects 

We now turn to a qualitative discussion of interaction ef- 
fects on the system. To do this we employ Bogoliubov's 
method for treating an interacting Bose gas. In this approach 
all the atoms are assumed to be in the zero quasi-momentum 
condensate, and by diagonalizing an approximation to the 
many-body Hamiltonian a set of quasiparticle levels are ob- 
tained (e.g. see 1121 lOO . This approach should be a good 
approximation to the exact many-body spectrum for moderate 
lattice depths (i.e. below the depth where the Mott-transition 
occurs) and small thermal depletion. 

The regimes of cooling associated with lattice loading ob- 
served in Figs. 12141 arose from the rapid compression of the 
ground band width (ebw) that occurs with increasing lattice 
depth. To assess the effects of interactions on cooling we com- 
pare the width of the ground band calculated with Bogoliubov 
theory to the non-interacting case for parameters relevant to 
current experiments in Fig. 0a). These results demonstrate 
that interactions suppress the rate at which the bandwidth de- 
creases as the lattice is applied. We also note that the band 
gap is modified by interactions, however the thermodynamic 
properties of the system are not as sensitive to small changes 
in this and we will not consider the effect of this. 
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Figure 5: Bogoliubov description of interaction effects on the (a) 
bandwidth and (b) scaling parameter a for various lattice depths. 
Crosses: interacting result for n = 1. Circles: non-interacting result. 
Other parameters: A = 850nm and collisional interaction strength 
taken forRb 87 . 



To quantify how the modified ground band excitation spec- 
trum affects thermodynamic properties of the system we con- 
sider an extension of the arguments made in Sec. IIII B 21 In 
that section we discussed the scaling relationship that could 
be applied to the thermodynamic quantities for a Bose gas in 
the tight-binding regime when only states of the ground band 
were relevant. For the interacting case we can apply the same 
argument if we ignore any ground band spectrum reshaping 
with varying lattice depth, other than a scale change given by 
the scaling parameter a — (eBw)//(eBw) t . where (eBw)i 
and (cbw)/ are the respective bandwidths of the initial and 
final configurations (both assumed to be in the tightbinding 
limit). We note that this definition of a is equivalent to Eq. 
H2i for the non-interacting case. As was shown in Eq. i ll Oi l, 
the ratio of final to initial temperatures is equal to a. In Fig. 
|5jb) we compare the values of a both with and without inter- 
actions for a system initially in a V = 4Er deep lattice, for 
various final depths V 11% . These results quantify how inter- 
actions lead to a lesser degree of cooling, and (in the limit of 
the Bogoliubov approach) these results suggest that the tem- 
perature will level out to some finite non-zero value at large 
lattice depths. 

In the Mott-insulator regime (where Bogoliubov approach 
is invalid) an energy gap proportional to the on-site interaction 
strength develops between the ground state and the particle- 
hole excitations above it (see 1 14]). This suggests that in gen- 
eral interactions tend to increase ground band width, and will 
limit our ability to cool the system with adiabatic passage. 



E. Adiabaticity 

Finally we note that interactions between particles are es- 
sential for establishing equilibrium in the system, and un- 



6 




Figure 6: Fast lattice loading of a N s ~ 3 x 10 4 site cubic lat- 
tice, with filling factor n = 1. Broken/dottedVdashed lines indicate 
fast loading curves (see text) and are labeled by the temperature of 
the initial V = state. The lattice depth on these curves can be 
determined from their intercept with the equilibrium entropy versus 
temperature curves (gray solid lines), which are described in Fig. [3] 



derstanding this in detail will be necessary to determine the 
timescale for adiabatic loading. In general this requirement 
is difficult to assess, and in systems where there is an addi- 
tional external potential it seems that the adiabaticity require- 
ments will likely be dominated by the process of atom trans- 
port within the lattice to keep the chemical potential uniform, 
though recent proposals have suggested ways of reducing this 
problem fl5ll . 

It seems reasonable that for sufficiently deep lattices the 
decreasing tunneling rate will ultimately become the rate lim- 
iting timescale for maintaining adiabaticity in adiabatic load- 
ing (e.g. see 0). To estimate this timescale we consider 
a case relevant to 87 Rb experiments. For a lattice depth of a 
V = lOEji, where we have taken A = 850nm, the tunneling 
time is ~ 10ms. This timescale is short compared to the load- 
ing time used in recent experiments with this system (jlHI]' 
and suggests that smoothly increasing the lattice to depths of 
~ IOEr over > 100ms should be very adiabatic. For depths 
larger that V = IOEr the tunneling time increases exponen- 
tially and the adiabatic condition becomes more difficult to 
satisfy, however this is also the regime where many-body ef- 
fects will begin to dominate and a more complete description 
will be needed to fully understand the adiabaticity require- 
ments. 

It is useful to assess the degree to which non-adiabatic load- 
ing would cause heating in the system. We consider lattice 
loading on a time scale fast compared to the typical colli- 
sion time between atoms, yet slow enough to be quantum me- 
chanically adiabatic with respect to the single particle states. 
This latter requirement excludes changing the lattice so fast 
that band excitations are induced, and it has been shown that 
in practice this condition can be satisfied on very short time 



scales fl7ll . We will refer to this type of loading as fast lattice 
loading, to distinguish it from the fully adiabatic loading we 
have been considering thus far. 

To simulate the fast lattice loading we take the system to 
be initially in equilibrium at temperature Ti for zero lattice 
depth. For the final lattice depth we fast load into, we map the 
initial single particle distribution onto their equivalent states 
in the final lattice, and calculate the total energy for this fi- 
nal non-equilibrium configuration. This procedure assumes 
that there has been no collisional redistribution to allow the 
system to adjust to the changing potential. To determine the 
thermodynamic state the final distribution will relax to, we use 
the energy of the non-equilibrium distribution as a constraint 
for finding the equilibrium values of temperature and entropy. 
In general the final state properties will depend on the initial 
temperature, filling factor, and final depth of the lattice, and to 
illustrate typical behavior we show a set of fast loading pro- 
cess curves that indicate the final state equilibrium properties 
as a function of final lattice depth in Fig. [6]for unit filling and 
various initial temperatures. 

These curves show, as is expected from standard thermo- 
dynamic arguments, that entropy always increases for non- 
adiabatic processes, i.e. all loading curves in Fig. [6] bend 
up wards with increasing lattice depth. We also observe that 
eventually all curves predict heating for large enough final lat- 
tice depth. However, for initial temperatures sufficiently far 
below the critical temperature for cooling (Tq ~ 0.7En/kB) 
a useful degree of temperature reduction can be achieved with 
fast lattice loading up to certain maximum depth. For exam- 
ple, the curve with initial temperature Ti = O.^Er/Ice in 
Fig. [6] cools for final depths less than V ~ 5Er, but for 
depths greater than this the temperature begins to increase 
quite rapidly. 

IV. 2D LATTICE WITH DIPOLE CONFINEMENT 

Here we consider the case of a 2D lattice in the xy-plane 
with a harmonic trapping potential along the z direction. We 
calculate the eigenstates of the Schrodinger equation with po- 
tential energy 

V T {v) = j [cos(2kx) + cos(2fcy)] + ^m(2ir f T ) 2 z 2 , 

where fx is the harmonic trap frequency. In Fig. 0we give the 
entropy versus temperature curves for a range of experimen- 
tally relevant parameters. The results in Figs. 0(a) and (c) 
demonstrate the behavior for a range of lattice depths and two 
different atom densities, but with constant harmonic trap fre- 
quency. These curves exhibit qualitative similarities with the 
3D lattice curves shown in Figs. 12 141 such as a crossing re- 
gion where the zero-depth curve changes from being a lower 
bound to an upper bound of the other curves, and separates 
the regions where adiabatic loading of the system will cool or 
heat the system respectively. Similar to the 3D case, this be- 
havior arises because of how the lattice modifies the energy 
spectrum, i.e. the compression of low lying states and upward 
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Figure 7: Entropy versus temperature curves for a N s ~ 10 2 site 2D 
square lattice with harmonic trapping potential in the perpendicular 
direction, for various lattice depths and strengths of harmonic con- 
finement. Results are given 2D site occupations of (a)-(b) n = 10, 
and (c)-(d) n — 100. In each subplot the various curves correspond 
to lattices depths of V = to 20Er (with a spacing of 4E between 
each curve). In (a) and (c) the harmonic confinement frequency is 
f T = Q.QlEn/h, whereas in (b) and (d) f T = 0.0lE R /h at V = 0, 
and is taken to increase by 2 x 10 _3 -Eh for each subsequent lattice 
depth, to /t = 0.02£'_r at V — 20Er. The condensation point for 
the system is marked with a circle on each curve. 



shift of high lying states. However, the availability of equally 
spaced harmonic oscillator states ensures that the density of 
states does not have a gap (assuming fx is small compared to 
V) and an entropy plateau does not appear. We note that the 
condensation temperature decreases more rapidly both with 
temperature and entropy compared to the 3D cases, making 
the 2D lattice a more ideal system for observing reversible 
condensation. 

In Figs. 0(b) and (d) we consider equivalent systems to 
those in Figs. 0(a) and (c), except that we increase the trap 



frequency with lattice depth to model the effects of additional 
dipole confinement on the system. For the case of rubidium 
and taking A ~ 805nm, these results correspond to a harmonic 
confinement of fx ~ 35Hz at zero lattice depth, with the con- 
finement increasing in linear steps on successive curves up to a 
maximum value fx ~ 70Hz at V — 2QEr, typical of current 
experimental parameters. These results demonstrate that the 
additional dipole confinement reduces the size of the region 
over which cooling occurs, and reduces the extend to which 
the system can be cooled. The more rapidly fx increases with 
lattice depth, the more pronounced this reduction will be. 



V. CONCLUSION 

In this paper we have calculated the entropy-temperature 
curves for bosons in a 3D optical lattice, and a 2D lattice 
with harmonic confinement for various depths and filling fac- 
tors. We have identified general features of the thermody- 
namic properties relevant to lattice loading, indicated regimes 
where adiabatically changing the lattice depth will cause heat- 
ing or cooling of the atomic sample, and have provided lim- 
iting results for the behavior of the entropy curves. We have 
considered the effect of lattice depth and filling factor on the 
Bose condensation point and have examined the possibility of 
reversible condensation through lattice loading. We have dis- 
cussed the dominant effects of interactions, and have shown 
that many of our predictions are robust to non-adiabatic ef- 
fects. Future extensions to this work will consider in more 
detail the effects of both interactions and inhomogeneous ex- 
ternal potentials. 
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